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Abstract 

We numerically investigate the area statistics of the outer boundary of planar 
random loops, on the square and triangular lattices. Our Monte Carlo simulations 
suggest that the underlying limit distribution is the Airy distribution, which was 
recently found to appear also as the area distribution in the model of self-avoiding 
loops. 

1 Introduction 

In recent years, a number of predictions for critical exponents related to random planar 
curves have been proved using probabilistic methods. Particular examples are the proof 
of old conjectures concerning exponents of critical percolation on the triangular lattice £Q, 
and the proof of the conjecture that the Hausdorff dimension of the outer boundary of 
planar Brownian motion is 4/3 2\. The latter conjecture was stated in 1982 on the basis of 
an investigation of random loops on the square lattice j3] , while a non-rigorous derivation 
using methods of quantum gravity was given in 1998 [I]. 

A main ingredient of the approach ^ |2] is a family SLE K of conformally invari- 
ant stochastic processes, called Schramm-Loewner evolution (sometimes also stochastic 
Loewner evolution). Critical percolation and the outer boundary of planar Brownian mo- 
tion are described by SLEq. More generally, it is conjectured that Schramm-Loewner 
evolution describes critical Fortuin-Kasteleyn clusters joj- The area distribution of such 
clusters has been analysed recently 0. 

Schramm-Loewner evolution is also connected to the self-avoiding walk model jjj. It 
has recently been shown that, if the continuum limit of self-avoiding walks exists and is 
conformally invariant, it is described by SLEg/3 [Bj. This led to another derivation of the 
critical exponents of the self-avoiding walk v — 3/4 and 7 = 43/32, and to predictions of 
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the distribution of some random variables related to the self-avoiding walk problem. Some 
of these have been numerically tested recently Moreover, there are predictions for 

the self- avoiding polygon model |J] based on SLEg/3. Essentially, the scaling limit of self- 
avoiding polygons is predicted to be described by the union of two Brownian excursions 
[Hj. Also, it is predicted that the scaling limit coincides with the outer boundary of loops 
of Brownian motion [H]. Since the area under a Brownian excursion is Airy distributed 
[TT] , one thus expects the Airy distribution 1 as the area distribution of the outer boundary 
of Brownian motion. 

A completely different, independent approach to the self-avoiding polygon problem is 
inspired by the analysis of exactly solvable models of lattice polygons, counted by perimeter 
and area. A conjecture for the scaling function of self-avoiding polygons [T3J[T3] has recently 
been tested [T3J E3 QUI, on the basis of an extrapolation of exact enumeration data and 
Monte Carlo simulations on the square, triangular and hexagonal lattices. The validity of 
the conjecture implies that the asymptotic distribution of area of self-avoiding polygons is 
the Airy distribution. 

Within this framework, the conjecture for the scaling function of self-avoiding polygons 
arose from an analysis of polygon models whose perimeter and area generating function 
satisfy a g-algebraic functional equation compare also ^7]. It was shown that if 
their perimeter generating function displays a square-root singularity as the dominant 
singularity, then the mean area a m of a polygon of perimeter m generically grows as a m ~ 
Aw? v as m — » 00, where v = 3/4. This implies a dimension \ jv = 4/3 for such polygons. 
Moreover, a certain scaling function of Airy- type generically arises fH] • Since the prediction 
of the dimension of self-avoiding polygons 4/3 was also supported by conformal field theory 
arguments and numerical analysis j7j , it was natural to ask whether self-avoiding polygons 
might have the same scaling function. Noting that the outer boundary of planar random 
loops also has dimension 4/3, one might suspect that the Airy distribution appears here 
as well. 

In this article, we test the conjecture that the Airy distribution arises as the limit dis- 
tribution of the outer boundary of planar random loops, based on Monte Carlo simulations 
on the square and triangular lattices. Our conjecture arises from the connection to the 
self-avoiding polygon model, which may be viewed from either the stochastic perspective 
or from the exactly solvable polygon model perspective, as indicated above. We find the 
conjecture of an Airy distribution for the limit distribution of the outer boundary of planar 
random loops satisfied, within numerical accuracy. 

The article is organised as follows. In the following section, we will gather basic results 
about the Airy distribution. The third section explains our algorithms and presents the 
analysis of data. We then summarise our results and discuss possible future work. In 
the appendix, we discuss the connection between scaling functions and limit distributions. 
This motivates the particular form of our prediction for the area moments in section two. 

1 The Airy distribution appears in a variety of other contexts such the path lengths in trees and the 
analysis of the linear probing hashing algorithm, see |12j for a review. 
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2 Airy distribution 



A random variable X is said to be Airy distributed if its moments are given by 

where the Airy constants flk are determined by the quadratic recurrence 

k—i /,\ 

O = -1, 20 fe = (3k - 4)Jfefi fc _i + ( • ) ^jfifc-i (*>!)■ (2) 

The first few numbers tt k are -1, 1/2, 5/4, 45/4, 3315/16, 25425/4, 18635625/64, 18592875. 
The moments uniquely define a probability distribution, which is called the Airy distri- 
bution. The name relates to the fact that the numbers Qk also appear in the asymptotic 
expansion 

^logAiM-fi^!-^ (3 ) 

A:=0 

where Ai(z) = - J °°cos(t 3 /3 + tz) is the Airy function. 

In the following section, we will analyse the conjecture that the outer boundary of planar 
random loops is Airy distributed, using a Monte Carlo simulation, both for the square and 
triangular lattices. More precisely, denote by p m ,n the number of sampled outer boundaries 
of planar random loops with perimeter m and area n. (An algorithmic definition of the 
outer boundary of a random loop will be presented in the following section.) Define a 
corresponding random variable Y m via 

HYm = n) = (4) 

2_^ n Pm,n 

We will numerically analyse the area moments E[Y^] under the assumption that 

for some constant C > 0. In comparison to ([Tjl. the factor m 3k ^ 2 accounts for the growth of 
the area moments in powers of m 3//2 . This formula correctly describes the area statistics of 
the self-avoiding polygon model on the square and triangular lattice, with p m ,n the number 
of (sampled) self-avoiding polygons of perimeter m and area n. (The perimeter of a polygon 
is the number of edges, its area is the number of enclosed squares or triangles.) See the 
appendix for a derivation of (J3J) from the conjectured scaling function of self-avoiding 
polygons. 

2 For the definition of the Airy distribution and its properties we follow |12| . 
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3 Numerical analysis 



Let us consider random walks on a planar lattice. For those walks returning to their 
starting point, we construct the outer boundary using the following three-step algorithm. 

(1) Find the bottom vertex Vq of the random walk. This is the first vertex in an increasing 
lexicographic ordering of the walk vertices by their coordinates; first in the ^-direction 
and then in the y-direction. Define the direction vector e Vo = —e y , where e y is the 
unit vector in the ^-direction. Set i = 0. 

(2) For the vertex V{, consider the set Fi of direction vectors pointing from Vi to those 
random walk vertices which are also lattice nearest neighbours of t>j. Choose from Fi 
the direction vector fi of maximal angle a(e Vi , fi) G [0, 2tt). Set t>j + i = Vi + fi and 

e v i+1 = ~fi- 

(3) If v i+ i 7^ vq, then repeat step (2), with i replaced by i + 1. Otherwise, set N = i and 
stop. 

The vertices V{, subsequently obtained from the algorithm, yield a closed walk (t> , v%, . . . , vn), 
called the outer boundary of the random walk. We define its perimeter to be N + 1. The 
walk travels along the random walk in clockwise direction. A typical square lattice example 
is shown in Fig. ^ Note that the outer boundary is not a self-avoiding polygon, since it has 




Figure 1: A square lattice random loop (left) and the corresponding outer boundary (right). 

cut vertices. The simulations suggest however that such vertices occur with zero density 
in the limit of infinite length of the outer boundary. 

We considered the cases of the square lattice and the triangular lattice. The triangular 
lattice case was implemented on the computer using a square lattice with additional edges 
along the positive diagonal of each square. Due to piecewise linearity, the area A of a loop 
(vq, v i, . . . , vn) is, on both lattices, given by 

N 

A = 2^y Vi + 1 > x ~ — 2 — ' 

i=0 
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where Vi jX and v^ y denote the x-coordinate and y-coordinate of the vertex u i( and vn + \ = Vq. 
It can thus be determined "on the fly" together with the construction of the outer boundary. 

We generated random walks up to length 10 5 on the square and triangular lattices, 
using the random number generator ran2 . c 19|. For each lattice, we generated 10 8 random 
walks. The size of the array was chosen such that the random walks did not reach the array 
boundaries. In our situation, an array of size 5000 x 5000 was sufficient. For those random 
walks returning to their start point, we constructed the outer boundary and computed the 
values of the first ten area moments for perimeter up to 10 3 , i.e., we sampled E[Y^] for 
m < 10 3 and k < 10. Note that only even values of m occur for the square lattice. For 
the data analysis, we rejected all values of m with less than 10 4 samples. Furthermore, we 
chose m > 50 in order to reduce finite size effects. This led to 50 < m < 512 for the square 
lattice and to 50 < m < 401 for the triangular lattice. The total number of loops used 
in the data analysis was 8.6 x 10 6 for the square lattice and 11.6 x 10 6 for the triangular 
lattice. 

We first analysed the predicted power-law behaviour of the area moments (jSJ) 

E[Y£] ~ A k m a " (m -> oo), (6) 

with cx,k — 3/c/2. Estimates of a k were obtained by a least squares fit of the data to the 
assumed asymptotic form 

N B 

log(E[F^]) ~ log A k + a fc logm + ^ -JL ( m _ oo). ( 7 ) 

Since the exponents of the corrections to the asymptotic behaviour of © are not known, 
we chose /3k,i = I- This corresponds to analytic correction terms in (j^J). The results of 
a least squares fit with N = 3 are shown in Table ^ The numbers in brackets give the 
statistical error of the fit in the last two digits. We found that differences in the estimates 
of a k due to different choices of the exponents ftk,i are an order in magnitude smaller than 
the statistical errors. This reflects the observation that the statistical fluctuations of the 
data points dominate effects due to different corrections to the asymptotic behaviour (JHJ). 
The obtained values are consistent with the prediction ak = 3fc/2, which we employ in the 
following. 

We then tested our prediction for the area moment amplitudes. It follows from (jSJ) that 
the ratios AkjA\ are independent of the constant C. In particular, we get 

A*/ A? = 2 — k [gjLgj;^*, A 2k+1 /Al k+1 = ^^W, (8) 

where the numbers Qk are the Airy constants defined in (J2J). We extrapolated the numbers 
Ak/A\ by a least squares fit to E[F^]/E[F m ] fc ~ A k /A\ + c k /m + d k /m 2 as m —> oo. A 
typical data plot is shown in Figure El 

The corresponding results for the first ten amplitude ratios are given in Tabled together 
with their predicted values. The numbers in brackets give the statistical error of the fit in 
the last two digits. 
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k 


1 


2 


3 


4 


5 


Square ak 
Triangular ak 


1.522(41) 
1.508(81) 


3.043(85) 
3.02(16) 


4.56(14) 
4.53(24) 


6.08(20) 
6.04(34) 


7.60(27) 
7.56(45) 



k 


6 


7 


8 


9 


10 


Square ak 


9.12(36) 


10.63(48) 


12.14(62) 


13.65(80) 


15.1(10) 


Triangular ak 


9.08(58) 


10.60(73) 


12.11(90) 


13.6(11) 


15.2(14) 



Table 1: Estimates of the exponents ak of the area moments E[Y^] on the square and 
triangular lattices. The numbers in brackets denote the statistical error of the fit in the 
last two digits. 



Amplitude 


Airy distribution 


Square MC 


Triangular MC 


A 2 /A\ 


1.06103296 


1.0626(18) 


1.0620(23) 


As/Al 


1.19366208 


1.1980(52) 


1.1951(61) 


A 4 /Af 


1.42171312 


1.430(13) 


1.422(15) 


As/Al 


1.78895216 


1.801(28) 


1.783(31) 


A % /A\ 


2.37205758 


2.387(59) 


2.350(61) 


A 7 /A\ 


3.30494239 


3.32(12) 


3.24(12) 


As/Al 


4.82446133 


4.83(25) 


4.68(24) 


A 9 /Al 


7.35752352 


7.32(52) 


7.02(48) 


A 10 /A{° 


11.6901282 


11.5(11) 


10.93(98) 



Table 2: Comparison of amplitude ratios arising from the Airy distribution against MC 
estimates for the square and triangular lattice. The numbers in brackets denote the sta- 
tistical error of the fit in the last two digits. 
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EK]/E[K 



12 1.12 



1.11 



1 . 1 



1 .09 



1 .08 



1 .07 




O 




O.Ol 0.02 0.03 0.04 



1/m 



Figure 2: MC estimates of E[Y*]/E[Y m } 2 on the square lattice, plotted against 1/m. 



Finally, we estimated the constant C in (jSJ) using the first area moment K[Y„ 
7rm 3 / 2 as m 



oo. A least squares fit to E[Y m ]/m 3 / 2 with a polynomial of degree 4 in 
1/m yields C sq = 0.0910(11) for the square lattice and C tr = 0.1262(29) for the triangular 
lattice. 



4 Conclusion 

Our numerical analysis supports the conjecture that the area of the outer boundary of 
planar random loops is Airy distributed, on both the square and triangular lattices. Since 
we expect this result to be universal, it would be interesting to check the appearance of the 
Airy distribution for boundaries of random walks on, for example, non-periodic tilings such 
as the Penrose tiling or the Ammann-Beenker tiling. Also, the question arises whether the 
Airy distribution can be found in other models of cluster boundaries with fractal dimension 
4/3. The outstanding question is whether our conjecture can be proved, possibly by using 
quantum gravity methods [I] or the probabilistic approach j2]. 
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Appendix: Scaling functions and limit distributions 

We explain how limit distributions are derived from the scaling function of a combinatorial 
model. We show, in particular, how the Airy distribution is derived from the scaling 
function of polygon models [2H] such as staircase polygons or self-avoiding polygons. This 
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relates the scaling function viewpoint employed in [T3J EH UH] to the present approach. In 
particular, it motivates formula (J5J). 

For a polygon model on a given lattice, denote the number of polygons of (half-) perime- 
ter m and area n by p m ,n- The perimeter and area generating function of the model is 
given by G(x,q) = J2 mn Pm.,nX m q n . The function G(x, 1) = go(x) is called the perimeter 
generating function. Let x c denote the radius of convergence of the perimeter generating 
function. Typically [18J, the factorial area moment generating functions gk{x) have the 
asymptotic behaviour 



l) k d 
dq 



k G{x,q) 



fk 



9=1 



(x c - x)"»< 



x„ 



(9) 



with critical amplitudes fk and critical exponents jk, where 7^+1 > -fk- See jSHl Ch. 7] for 
a proof that the functions gk(x) exist and are finite for polygon models. It is convenient 
to incorporate the numbers fk and 7*, into the function (THl HE] 



F(s) 



00 „ 

E Jk 



k=0 



(10) 



Typically, the function F(s) describes scaling behaviour of the perimeter and area gener- 
ating function about (x, q) = (x c , 1), 



G (sm9) {x,q) ~ (l-q) 9 F 



[1-qY 



(x,q) 



[x c , 1" 



'11} 



with critical exponents 9 and 0. The function F(s) is also called the scaling function. The 
form (jlljl is consistent with © if the exponents 7^ are of the form 7^ = (k — &)/<fi- One 
should note that (jlip has been proved only for a few examples including staircase polygons 
[2*Tj . For typical polygon models such as staircase polygons 3 and rooted self-avoiding 
polygons PHI, the scaling function has the particular form ^SJ Eqn. 3.9] 



F(s) = -Af 1 — logAi 
ds 




with critical exponents 9 = 1/3 and <ft — 2/3. 

We now turn to a probabilistic description, see also Selecting polygons uniformly 
among the set of all polygons with perimeter m defines a discrete random variable X m 
with probability distribution 



n 



m,n 



(13) 



3 More generally, for polygon models satisfying a q-algebraic functional equation with a square-root 
singularity as dominant singularity of their perimeter generating function |18| . 
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whose moments are given by 

ELY*] = ^ nkpm >\ (14) 

Z^n Pm,n 

The numerator in (JTljl is related to the coefficients of the factorial moment generating 
functions © via 

n n 

(15) 

where (n)^ = n(n — 1) • • • (n — + 1), and T(z) is the Gamma function. This relation 
can be used to characterise the leading asymptotic behaviour of the moments ELY^J of 
the underlying distribution. Let us introduce distribution coefficients defined by fk = 



^2 fc $ fe / 1 fc / 1 - fe . We obtain 



V0 x k 



^-M) (I) f * (16 > 



Now introduce a normalised random variable Y m defined by 

= "w- (17) 

„ Xc J fo 

The moments E[Y^J of the normalised random variable X m are then asymptotically given 
by 

nx k j ~ -^f* ("» - «>)■ (is) 

Note that (JTSj) is independent of x c and asymptotically constant. In conjunction with Q, 
it follows that the particular scaling function f!12D leads to the Airy distribution (0) in (|18j) . 
Formula © is chosen in analogy with pfij) . 
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